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Two methods are described to obtain long-term (>20 years) distri- 
butions of 5-minute point rain rates from data published by the Na- 
tional Climatic Center for U.S. locations. A set of simple empirical 
formulas for converting the distribution of 5-minute rain rates into rain 
attenuation distributions on 11 -GHz radio paths has been deduced 
from data measured in Georgia. Additional data measured in several 
other locations also support this empirical formulation. These simple 
formulas and 5-minute point rain rate distributions are useful for path 
engineering of 11 -GHz radio. The work on rain rate distributions dis- 
cussed in the paper derives from approaches suggested by the late W. 
Y. S. Chen. 5 - 9 ' 16 

I. INTRODUCTION 

An important problem in designing terrestrial and earth-satellite radio 
systems at frequencies above 10 GHz is the added path attenuation 
caused by rain. Long term (>20 years) rain rate statistics are needed to 
engineer radio paths for various geographic locations to meet reliability 
objectives. Section II describes a method to obtain 20-year distributions 
of 5-minute point rain rates from the excessive short duration rainfall 
data 1 - 2 for locations in the relatively wet eastern and midwestern U.S.A. 
Section III describes another method, employing the theory of extreme 
value statistics, to obtain 50-year distributions of 5-minute point rain 
rates from rainfall-intensity-duration-frequency curves 3 for locations 
in relatively dry locations such as western U.S.A. Section IV discusses 
the variability of rain rate distributions with observation time base. 
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Table I — Thresholds* of excessive short-duration rainfalls 





Minimum 


Threshold 






depth of 


(r minute 






recorded 


average 




Duration t, 


rainfall, 


rain rate), 




minutes 


inches 


mm/hr 




5 


0.25 


76.2 




10 


0.30 


45.7 




15 


0.35 


35.6 




20 


0.40 


30.5 




30 


0.50 


25.4 




45 


0.65 


22.0 




60 


0.80 


20.3 




80 


1.00 


19.1 




100 


1.20 


18.3 




120 


1.40 


17.8 




150 


1.70 


17.3 




180 


2.00 


16.9 





* Exceeding any one of these 12 thresholds is sufficient to qualify a rainstorm as an 
excessive rainfall. Therefore, this definition does not require that an excessive rainfall 
exceeds all the 12 thresholds. 

In this paper, a "5-minute rain rate" corresponds to the average value 
of the randomly varying rain rate in a 5-minute interval and is calculated 
as AH/t where AH is the 5-minute accumulated depth of rainfall and 
t = 5 minutes = V12 hour is the rain gauge integration time. Similarly, 
a "r-minute rain rate" is the average rain rate in a r-minute interval. 

Sections V to VII describe a set of simple empirical formulas deduced 
from experimental data in Georgia for converting the distribution of 
5-minute point rain rates into distributions of rain attenuation on 11- 
GHz radio paths. Section VIII compares the calculated results with 
additional data from other locations. 

II. 20-YEAR DISTRIBUTIONS OF 5-MINUTE RAIN RATES 

The excessive short duration rainfall data 1 record details of those 
heavy rainfalls which exceed one or more thresholds; these thresholds 
are dependent upon the rain integration times t as shown in Table I. For 
example, the thresholds are 76 and 20 mm/hr for r = 5 and 60 minutes, 
respectively. The data, published in tabular form, consist of a storm- 
by-storm compilation of accumulated depth of rainfall in the most in- 
tense 5, 10, 15, 20, 30, 45, 60, 80, 100, 150, and 180 minute periods. For 
example, Table II shows such data for Newark, New Jersey in 1972. Only 
three rainstorms exceeded the excessive rainfall thresholds at Newark 
during 1972. More detailed discussions on the excessive short duration 
rainfall data can be found in Refs. 1, 5 and 6. 

The method for obtaining 5-minute rain rates from this data source 
is illustrated in Table III for the storm of August 26, 1972 at Newark. In 
essence, for each storm, the most intense 5-minute accumulation gives 
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one sample of 5-minute rain rate; the difference between the 10-minute 
and 5-minute greatest accumulations gives the second sample of 5- 
minute rain rate; and so on. By applying this single operation to 20 years 
(1953 to 1972) of such data for Newark, New Jersey, we obtain a 20-year 
distribution of 5-minute rain rates in the range above the 76 mm/hr 
threshold as shown in Fig. 1. The results of this method were tested 5 
using the reports of the 20 storms which were available in more detailed 
form 6 and excellent agreement was found. 5 



10' t= 



< 
> 
IE < 

££ 

^< 

S hi 

feo 

u X 

cc uj 

^< 

z" 10° 



10" 



NEWARK, N.J. 
20 YEARS (1953-1972) 



J L 



J I I L 



20 40 60 



80 100 120 140 160 180 200 220 240 260 280 300 
5-MINUTE RAIN RATE, mm/hr 



Fig. 1 — Twenty-year distribution of 5-minute rain rates above the 76-mm/hr threshold 
at Newark, New Jersey. 



As observed in Table I, the published 5-minute rain rates are given 
only for those storms which exceed the thresholds in Table I. These data 
are, therefore, incomplete with respect to 5-minute rain rates of less than 
76 mm/hr. However, in the same publication, progressively lower 
thresholds (see Table I) are used for longer durations; for example, the 
threshold at 60-minute duration is 20 mm/hr. The method developed 
to extend the distributions of 5-minute rain rate to very low rain rate 
employs measured distributions of 1-hour rain rates as the basic data, 
since these are readily available. 1>2 ' 7 By using long-term (>10 years) data 
from New York City, Miami, and McGill Observatory in Canada, a 
simple empirical formula was deduced for converting the 1-hour rain 
rate distribution into the 5-minute rain rate distribution in the low rain 
rate region. 8 A normalizing procedure is used in this conversion formula 



1584 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1977 



to account for the difference in rain characteristics between geographic 
areas. Twenty-year distributions of 5-minute rain rates for locations in 
the eastern and midwestern U.S.A. have been obtained by this process. 
Some results are given in Refs. 8 and 9. 

In low rain rate areas such as Oregon and Washington, however, few 
rainfalls exceed the criterial thresholds, and hence few are included in 
the excessive short duration rainfall data. For example, at Spokane, 
Washington, the 5-minute rain rate exceeded the 76 mm/hr threshold 
only once in the 20-year interval from 1953 to 1972. In such areas, a 
different data source and a method using statistics of extremes, discussed 
in the next section, is more suitable. 

III. 50-YEAR RAIN RATE DISTRIBUTIONS AND EXTREME VALUE 
STATISTICS 

The statistical behavior of the extremes of a random variable has been 
extensively investigated. 10-15 In an unpublished work, W. Y. S. Chen and 
R. L. Lahlum 16 have applied the theoretical distribution of yearly 
maximum 5-minute rain rates and an empirical extrapolation to obtain 
the rain rate distribution in the range of interest to radio path engi- 
neering. We extend Chen and Lahlum 's method by incorporating the 
theoretical distributions of the yearly K largest 5-minute rain rates for 
K ranging from 1 to 12. The application of the higher-order statistics of 
extremes eliminates the need for empirical extrapolation. 

Briefly, the average time per year that a rain rate r, measured by a 
gauge with integration time t, is exceeded, is given by 9 



T{R>r)^ T £ 



l-e-e~ y £ 



K-lp-Ny 



(1) 



N=0 N\ 

for the range T(R > r) < 50 minutes/year, where 

y = a[(ln r) - U] (2) 

is called the reduced variate, and a and U are scale and location pa- 
rameters, respectively. Notice that T(R > r) in eq. (1) is uniquely de- 
termined by the two parameters a and U. These two parameters can be 
calculated from the rainfall-intensity-duration-frequency curves 3 which 
are available for U.S. locations. These rainfall-intensity-duration-fre- 
quency curves are derived by the Gumbel method 11 - 12 using the theory 
of extreme value statistics and are based on approximately 50 years (1900 
to 1950) of rainfall data. 

From this data source, we need only the following three numbers for 
a given location to calculate the long-term distribution of 5-minute rain 
rates: 
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M = the number of years of rainfall data from which rainfall-inten- 
sity-duration-frequency curves are derived, 

r a = the extreme rain rate with 2-year return period, i.e., the rain rate 
which is exceeded once in 2 years, on average, by the yearly 
maximum 5-minute rain rates, 

fb = the extreme rain rate with 10-year return period, i.e., the rain 
rate which is exceeded once in 10 years, on average, by the yearly 
maximum 5-minute rain rates. 

The formulas for calculating a and U are: 9 

V6 
a = a a a z (3) 

IT 



Z jrl 

/ a/6 J 



where 



LT=[/„+^-[ 7 -- --%| (4) 

A a -A 



In r a - In n 



and 



<Tz 



(5) 



_ A a In r b - A b In r a 
A a - Ab 

K . _ ln U 10year8 ) . 0.3665 (7) 

\ 2 years — 1 year/ 

* — fafc. 10yBa ? ) »Mi (8) 

\ 10 years — 1 year/ 

7 = Euler's Constant « 0.5772 (9) 

zw= - ,n (- ln ]i7tr) (10) 

_ 1 M 

^iizo) (ii) 

_ 1 M 

^ 2 = T} E [Z(/)P (12) 



= y/Z*-Z 2 (13) 

For example, Fig. 2 shows a portion* of the rainfall-intensity-dura- 
tion-frequency curves for New York City. The required three numbers 
read from Fig. 2 are 



* The source curves in Ref. 3 cover a wider range for duration t from 5 to 1440 minutes 
and return period from 2 to 100 years. 
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20 30 40 50 100 

DURATION T, MINUTES 

Fig. 2— Rainfall intensity-duration-frequency curve by method of extreme values (after 
Gumbel) for New York City based on 49 years (1903-1951) of rainfall data. 

M = 49 years (1903 to 1951) 

r a = 4.4 inches/hr = 111.8 mm/hr 

Tb = 6.5 inches/hr = 165 mm/hr 

By substituting these three values into eqs. (3) to (13), we obtain 

a = 4.363 

U = 4.63 

Substituting this a, U pair into eqs. (1) and (2) yields the 49-year dis- 
tribution (dashed line) of 5-minute rain rates as shown in Fig. 3. The 
49-year (1903 to 1951) distribution obtained agrees well with the 20-year 
(1953 to 1972) distribution obtained by the method based on excessive 
short-duration rainfall data. 

Long-term distributions of 5-minute rain rates for U.S. locations can 
therefore easily be obtained by the extreme value method. 

IV. VARIABILITY OF SHORT-TERM DISTRIBUTIONS OF 5-MINUTE RAIN 
RATES 

Figure 4 shows that the 20-year distributions of 5-minute rain rates 
at Central Park,* La Guardia Airport,* and Newark Airport* agree 



* All within the New York Metropolitan area. 
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Fig. 3 — Comparison of 49-year (1903-1951) distribution of 5-minute rain rates calculated 
by extreme statistics method with measured 20-year (1953-1972) data at La Guardia 
Airport, New York City. 

closely. On the other hand, Figs. 5 and 6 show that 4-year (1969 to 1972) 
and 1-year (1972) distributions at these three locations differed signif- 
icantly. Figure 7 displays the convergence of the rain rate distributions 
at Newark as the time base is increased from 1 to 20 years. Figure 8 in- 
dicates that with a single rain gauge measurement, even a 20-year time 
base may still be insufficient to provide stable statistics of extremely high 
rain rates (beyond 160 mm/hr). 

Table IV lists the intervals, in minutes per year, that the rain rate 
exceeded 140 mm/hr at Newark, La Guardia Airport, and Central Park 
during the 20-year observation period. The last column contains the 
three-site summation for each year. The last three rows in Table IV in- 
dicate the 20-year average t, the difference At, between the worst year 
and the best year, and the ratio At ft respectively. Figure 9 shows that 
the normalized range of variations, At ft, increases rapidly as the rain 
rate increases from 80 to 160 mm/hr. This behavior is consistent with 
the divergence between the upper and lower envelopes in Fig. 7. In Fig. 
9, notice that the normalized range of variation, At ft, of individual sites 
are significantly greater than that of three-site summation for high rain 
rates. Since point rain rate statistics may be representative of a short 
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Fig. 4 — Twenty-year (1953-1972) distributions of 5-minute rain rates at three locations: 
La Guardia Airport, Central Park, and Newark Airport, in the New York Metropolitan 



radio hop, but three-site summation statistics are more representative 
for a multihop radio route, then, as shown in Fig. 9, the normalized range 
of variations, At/t, of annual outage time of a radio route can be expected 
to be much smaller than that of a short radio hop. In other words, the 
annual outage time of a radio route is statistically more stable than that 
of a short radio hop. Intuitively, the statistical stability of the accumu- 
lated outage time of a radio route stems from the partial compensation 
effect of the incoherent, random variations of individual-hop outage 
times as displayed in Figs. 5 and 6. 

The three-site, 20-year measurements yield 60 samples of annual 
accumulated time that rain rates exceed 140 mm/hr as listed in Table 
IV. The average value of these 60 samples is 0.97 minutes per year. Notice 
that 51 out of these 60 samples are less than the average value. Such 
nonsymmetric deviations of small sample data from long-term, large- 
sample average have already been observed and discussed in Ref. 17. 

These data indicate that rain rate statistics gathered from a single rain 
gauge measurement require a very long time base to yield stable statistics 
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Fig. 5 — Four-year (1969-1972) distributions of 5-minute rain rates at three locations: 
La Guard ia Airport, Central Park, and Newark Airport, in the New York Metropolitan 
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Fig. 9 — The effect of three-site summation on the normalized range of variations of 
annual distributions of 5-minute rain rate. At is the difference between the worst year and 
the best year in number of minutes per year rain rate exceeds abscissa during the 20-year 
period (1953-1972), and t is the number of minutes per year rain rate exceeds abscissa 
averaged over the 20-year period (see Table IV). 

Table IV — Number of minutes per year rain rate exceeds 

140 mm/hr 



Time base 



Newark 



La Guard ia 
Airport 



Central 
Park 



Three- 
site 
summation 



1953 


0.0 


0.0 


0.0 


0.0 


1954 


0.0 


0.0 


0.0 


0.0 


1955 


0.0 


0.0 


0.0 


0.0 


1956 


0.0 


0.0 


0.0 


0.0 


1957 


0.0 


0.0 


0.0 


0.0 


1958 


0.0 


0.0 


0.0 


0.0 


1959 


0.0 


0.0 


0.0 


0.0 


1960 


0.0 


0.0 


5.0 


5.0 


1961 


5.0 


0.0 


0.0 


5.0 


1962 


0.0 


0.0 


0.0 


0.0 


1963 


0.0 


0.0 


0.0 


0.0 


1964 


0.0 


0.0 


0.0 


0.0 


1965 


5.0 


0.0 


0.0 


5.0 


1966 


0.0 


5.0 


0.0 


5.0 


1967 


0.0 


0.0 


0.0 


0.0 


1968 


0.0 


0.0 


5.0 


5.0 


1969 


0.0 


13.3 


0.0 


13.3 


1970 


5.0 


0.0 


0.0 


5.0 


1971 


0.0 


6.7 


0.0 


6.7 


1972 


8.3 


0.0 


0.0 


8.3 


t (20-year 










average) 


1.2 


1.3 


0.5 


2.9 


At (worst 










year-best 


8.3 


13.3 


5.0 


13.3 


year) 

At 

7 










6.9 


10.2 


10.0 


4.6 
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Fig. 10 — Distributions of 11 -GHz rain attenuation measured on three paths near Atlanta, 
Georgia. 

for engineering a radio route. On the other hand, if the time base is not 
sufficiently long, the short term results tend to underestimate the 
long-term, large-sample average. 

V. ANALYSES OF MICROWAVE RAIN ATTENUATION 

For radio path engineering applications, a procedure is needed to 
calculate the rain attenuation distributions on microwave radio paths 
from the available rain rate distributions. Several independent theo- 
retical analyses relating rain attenuation distributions to radio path 
length have been developed. 18 - 2330 One analysis is based upon the ap- 
proximate log-normality of long-term distributions of rain attenuation 
and of point rain rates. 17-22 These two distributions are related by 
suitably derived parameters. Since existing theory for converting rain 
rate into rain attenuation applies to spatially inform rain rates, whereas 
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actual rainfalls are almost never uniform over a radio path, the hop is 
divided into incremental volumes in each of which the rain rate is uni- 
form; the total attenuation is then obtained by integrating over the path. 
Since the rain rates at various positions along the radio path only par- 
tially correlated, the increase of attenuation with radio path length is 
nonlinear. 
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Fig. 11 — Nonlinear dependence of 11 -GHz rain attenuation on path length measured 
at Palmetto, Georgia. 

In another analysis, 23 rain cells of circular cross section are assumed, 
allowing calculation of the probable length of path on which rain will fall 
from the probability of rain occurring at a point. It is found that the in- 
crease of attenuation with hop length is nonlinear because of the finite 
diameter of rain cells. In the limit, for hop lengths smaller than the cell 
diameter, the attenuation is almost proportional to path length, but for 
hops much larger than the cell diameter, it is almost independent of path 
length. 

Both analyses indicate that the nonlinearity of the path length de- 
pendence is a function of rain rate and rigorous derivations are fairly 
complex. For path engineering applications, it is desirable to have a 
simple empirical formula to describe this nonlinear behavior. The fol- 
lowing two sections describe the empirical formulas deduced from the 
available 11 -GHz rain attenuation data. 
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Fig. 12 — Rain rate distributions measured at Atlanta and Palmetto, Georgia. 



VI. RAIN ATTENUATION AND RAIN RATE DATA AT ATLANTA 

The 11 -GHz rain attenuation distributions obtained by simultaneous 
measurement on three paths (5.07, 25.6, and 42.5 km) near Atlanta 24 
from November, 1973 to July, 1974 are shown for illustration in Fig. 10. 
The path-length dependence derived from these data is indicated in Fig. 
11 which is a cross-plot relating attenuation observed to path length 
traversed, for fixed levels of probability. These data demonstrate that, 
for probability levels of 0.05 and 0.1 percent, rain attenuation increases 
nonlinearly with increased path length. 
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Fig. 13 — Nonlinear dependence of 11 -GHz rain attenuation on path length mea8ured 
at another location. 



Rain rate data for the above measuring interval are available from a 
tipping-bucket rain gauge at a common path terminal (Palmetto, 
Georgia) and from a Weather Bureau rain gauge at Atlanta Airport. 1,2 > 3 
Figure 12 shows the distributions of 5-minute point rain rates obtained 
from these two rain gauges, as well as long-term results from the same 
Atlanta station. 

Figure 13 shows another example of nonlinear increase of rain atten- 
uation with path length measured at another location. The data in Figs. 
13 and 15 are from the same city. 



VII. EMPIRICAL PATH-LENGTH DEPENDENCE 

Many authors 26-35 have pointed out that the relationship between the 
rain attenuation gradient, /3 in dB/km, and the point rain rate, R in 
mm/hr, can be approximately described by 

MR) = pR" (14) 

where the coefficient p and the exponent 77 depend on the radio fre- 
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Fig. 14 — Comparison of calculated (dashed lines) versus measured 11-GHz rain atten- 
uation distributions on 5.07, 25.6, and 42.5 km paths at Palmetto, Georgia. The calculations 
are based on 20-year distribution of 5-minute point rain rates at Atlanta, Georgia. 

quency, the polarization of the radio signal, and the canting angles of 
the oblate raindrops. Based on T. S. Chu's theoretical calculation 36 and 
the experimental data on the polarization dependence of rain attenuation 
in Refs. 37 and 38, the empirical formula for 11-GHz rain attenuation 
gradient is 



0v(R) = 0.0153/? 11909 dB/km 
for vertically polarized signals and 

Ph(R) = 0.0170A 1 - 2012 dB/km 



(15) 



(16) 



for horizontally polarized signals. 

If the rain rates were uniform over a radio path of length L (km), the 
path rain attenuation n(R,L) would be simply (3(R)L, but since actual 
rainfalls are not uniform, the increase of n(R,L) with L is nonlinear. In 
Ref. 39, a two-parameter empirical formula was shown to describe this 
nonlinear behavior. A single parameter variation also provides satis- 
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Fig. 15 — Comparison of calculated (dashed lines) versus measured 11-GHz rain atten- 
uation distributions on 11.3, 20.8, and 32.1 km paths at City A. The calculations are based 
on 20-year distribution of 5-minute point rain rates at City A. 



factory agreement: 



n(R,L) = 0(R)L 



(17) 



l+=r 



L(R) 



where the factor 



1 + =• 



(18) 



L(R) 

accounts for the partially correlated rain rate variations along the path, 
and L(R) is a characteristic path length such that the nonlinear factor 
(18) equals one-half when L = L. L is related to the diameter of the rain 
cell. 
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Fig. 16 — Comparison of calculated (dashed line) versus measured 11-GHz rain atten- 
uation distribution on a 10.1 -km path at City B. The calculation is based on 20-year dis- 
tribution of 5-minute point rain rate at City B. 



By fitting eq. (17) to concurrently measured (August 1973 to July 1974) 
distributions of 5-minute point rain rates and 11-GHz rain attenuation 
on the 42.5 km path at Palmetto, Georgia, it is found that L(R) can be 
approximately described by 



L(R) 



2636 
R-6.2 



km 



for 



R > 10 mm/hr 



(19) 



VIII. COMPARISON OF PREDICTED AND MEASURED 11-GHZ 
ATTENUATIONS 

Measured rain attenuation data on nine 11-GHz paths, listed in Table 
V, are available for comparison with the calculated results. Figures 14 
to 18 show such comparisons. The calculated results (dashed lines) are 
based on 20-year distributions of 5-minute point rain rates and include 
path rain attenuation and assumed wet radome loss listed in Table V. 
It is assumed that, on the average, a wet, flat radome introduces 2 dB 
loss and a wet, hemispheric radomes introduces 4-dB loss. 18 - 40 - 41 Figures 
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Fig. 17 — Comparison of calculated (dashed line) versus measured 11-GHz rain atten- 
uation distribution on a 12.2-km path at City C. The calculation is based on 20-year dis- 
tribution of 5-minute point rain rates at City C. 
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Fig. 18 — Comparison of calculated (dashed line) versus measured 11-GHz rain atten- 
uation distribution on a 32.3-km path at City D. The calculation is based on 20-year dis- 
tribution of 5-minute point rain rates at City D. 
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14 to 18 indicate that the measured results are comparable with the 
predicted curves. 

IX. CONCLUSION 

Two methods have been described to obtain long term (>20 years) 
distributions of 5-minute rain rates from data published by the National 
Climatic Center for U.S. locations. Some typical results are given in Refs. 
8 and 9. The variability of distributions based on shorter terms is dis- 
cussed in Section IV. 

A set of simple empirical formulas, for converting the distributions 
of 5-minute rain rates into the distributions of 11 -GHz rain attenuation 
on any path length, has been deduced from experimental data gathered 
near Atlanta, Georgia. These formulas are supported by further exper- 
imental data from other locations. 
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